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The PMNS and CKM matrices are phenomenologically close to symmetric, and a symmetric form 
could be used as zeroth-order approximation for both matrices. We study the possible theoretical 
origin of this feature in flavor symmetry models. We identify necessary geometric properties of 
discrete flavor symmetry groups that can lead to symmetric mixing matrices. Those properties are 
actually very common in discrete groups such as A4, S4 or A(96). As an application of our theorem, 
we generate a symmetric lepton mixing scheme with #12 = #23 = 36.21°; @13 = 12.20° and 5 = 0, 
realized with the group A(96). 


I. INTRODUCTION 


The properties of the fermion mixing matrices are ex¬ 
pected to give important hints on the underlying flavor 
physics. Flavor symmetries [T| are an attractive and 
most often studied approach to explain the rather differ¬ 
ent structure of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) and Cabibbo-Kobayashi-Maskawa (CKM) mix¬ 
ing matrices. Literally hundreds of models have been 
proposed in the literature, applying many possible dis¬ 
crete groups in order to explain lepton and quark mix¬ 
ing. Instead of adding simply another model to that list, 
we study in this paper an interesting possible property 
of both the CKM and PMNS matrices. Namely, despite 
the fact that the CKM mixing is a small while the PMNS 
mixing is large, both can to reasonable precision be es¬ 
timated to be symmetric. The symmetric form of the 
CKM matrix has early been noticed and studied in many 
references EHE]- After neutrino oscillation was well es¬ 
tablished the possible symmetric PMNS matrix also at¬ 
tracted some attention HME]. The symmetric form dis¬ 
cussed in these references includes the manifestly sym¬ 
metric case (U = U T ) and the Hermitian case (U = W). 
It is easy to get the relation 

(U = U T ) => {\U\ = \U\ T ) <=(!/ = U ] ) (1) 


by taking absolute values, which implies any physical pre¬ 
diction from \U\ = \U\ T can also be used in the other two 
cases U = U T or U = U'. Both of them are special cases 
of \U\ = \U\ T , which is what we mean by symmetric 
mixing matrix from now on. 

Using the global fits of the CKM [TS] and PMNS JT9] 
matrices, one finds: 


\U C 


( ( 0.97441 \ 
l 0.97413 ) 
( 0.22583 \ 
l 0.22461 ) 
( 0.00919 \ 
\ ^ 0.00854 ) 


( 0.22597 \ 
{ 0.22475 J 
f 0.97358 \ 
^ 0.97328 J 
0.0416 \ 
0.0393 J 


( 0.00370 \ 
y 0.00340 ) 
( 0.0426 \ 
^ 0.0402 ) 
( 0.99919 \ 
^ 0.99909 ) 
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/ 
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PMNS | 


/ 0.845 \ 
\ 0.791 J 
f 0.521 \ 
l 0.254 J 
( 0.521 \ 
^ 0.254 ) 


( 0.592 \ 
l 0.512 ) 
( 0.698 \ 
[ 0.455 J 
f 0.698 \ 
^ 0.455 J 


( 0.172 \ \ 
l 0.133 J 
f 0.782 N 
l 0.604 J 
f 0.782 \ 

^ 0.604 J ) 


( 3 ) 


Here the upper (lower) values in each entry are upper 
(lower) bounds of the matrix elements. The CKM matrix 
has been measured to a high precision (here we show the 
lcr range) and the relations \U\ 2 \ = IC/ 21 1, IC/ 23 1 = | U 32 \ 
are still well compatible with data. The relation | L/ 13 1 = 
|C /31 1 is, however, not fulfilled by data. As a symmetric 
mixing matrix requires that mm 


IC/ 31 I 2 - \U 13 \ 2 = | C / 12 | 2 - \u 21 \ 2 = \u 23 \ 2 - \u 32 \ 2 = 0 , 

( 4 ) 

we have an interesting option, namely, that some flavor 
symmetry or other mechanism generates \U\ 2 \ = |C/ 21 1, 

| C/ 23 1 = |C/ 32 1 but C /13 = C /31 = 0. Higher order correc¬ 
tions, which are frequently responsible for the smallest 
mixing angles, are then the source of non-zero | C/ 13 1 7 ^ 
|C/ 31 1, as well as for CP violation. Rather trivially, matri¬ 
ces with only one mixing angle are symmetric, the same 
holds for the unit matrix. 

The symmetry conjecture for the PMNS mixing is less 
compatible with data, as shown by the 3er bounds in Eq. 
© m- Similar to the quark sector, the 13- and 31- 
elements are incompatible with symmetry (the other two 
relations between the elements are also not favored by 
data), and a similar situation as mentioned above for the 
CKM matrix might be realized. Of course, one could also 
imagine that an originally symmetric mixing matrix is 
modified by higher order corrections, VEV-misalignment, 
RG-effects or other mechanisms that have been studied 
in the literature. 

For completeness, we give the phenomenological pre¬ 
diction of a symmetric mixing matrix, using the standard 
parametrization of the CKM and PMNS mixing matrices 

HJ: 


\U 33 \ = 


sin 0 i 2 sin 023 


v /1 — sin 2 S cos 2 0 i 9 cos 2 0 23 + cos 5 cos d\ 2 cos d 23 

( 5 ) 
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This is the unique physical prediction of both \U\ = \U\ T 
and U = U T . Note that \U\ = \U\ T has only one predic¬ 
tion as the unitarity requires that the relation in Eq. © 
is fulfilled, so once we set | C/ 13 1 2 = | C/ 31 1 2 we immediately 
get | U | = \U\ T . It is also the unique prediction of U = U T 
because any 3-by-3 unitary U with \U\ = \U\ T can be 
transformed to a new U' satisfying U' = U ,T simply by 
rephasing mm- Note that this does not hold for more 
than 3 generations. The Hermitian case U = U' has not 
only the prediction of Eq. (|5| but also CP conservation 
|13j . Thus, it would predict sin0 13 = tan0 12 tan 023- 

Despite the apparent deviation from their symmetric 
forms, one can still use it as an attractive zeroth-order 
Ansatz and attempt to study its theoretical origin. One 
option, put forward in m, is that a single unitary matrix 
V diagonalizes all mass matrices of quarks and leptons 
at leading order, in addition to the SU(5) relation rrid = 
mj between the down quark and charged lepton mass 
matrices. While it is difficult to embed this in realistic 
mass spectra in GUTs, the predictions of this scenario are 
that U G km = V^V = 1 while U PMNS = V T V is symmetric. 

In this paper, we study the the origin of symmetric 
mixing matrices from an underlying flavor symmetry. We 
prove a theorem which links geometric properties of dis¬ 
crete symmetry groups to the symmetric form of the mix¬ 
ing matrices. This theorem, explained in detail in Section 
PH holds only for subgroups of SO( 3) with real represen¬ 
tations, and can interestingly be realized in the most of¬ 
ten studied groups A 4 and S 4 . We find a modification 
that holds in the complex case in Section |III[ that could 
be applied to subgroups of 517(3), for instance to A(96). 
We use this to reproduce a previously studied, and actu¬ 
ally symmetric, mixing scenario for the PMNS matrix in 
Section EV] 

Since our analysis links the properties of the symmetry 
group with the mixing matrix, we end our introduction 
with a summary on how the generators of the group can 
be related to the matrices diagonalizing the mass ma¬ 
trices, following the strategy developed in In 

general, if a flavor symmetry group G is applied to, for 
instance, the lepton sector, then it must be broken to two 
residual symmetries Ge and G„ acting on the charged lep¬ 
ton sector and neutrino sector: 


G -a 


G v : 
G e : 


S T M V S = M v 
T^M(T = M t . 


( 6 ) 


Here the left-handed neutrino (assumed to be Majorana) 
mass matrix M„ is invariant under the transformation 
S T M V S for 5 £ G v and Me (defined by mem\, where me 
is the charged lepton mass matrix) is the effective mass 
matrix of left-handed charged leptons, invariant under 
T^MeT . Then the diagonalizing matrices U v and Ue de¬ 
fined by 

M v = U V D V U ?! (7) 

Me = U e DeUl ( 8 ) 

can be directly determined by 5 and T according to E3- 


U$SU„ = D s \ 
u\TUe = D t , 


( 9 ) 

( 10 ) 


where all D are diagonal matrices. Note that t/„ obtained 
from Eq. (JdJ) does not include Majorana phases which 
rephase each column of U u . Eq. 0 is independent of 
such rephasing which means the Majorana phases are 
not determined by flavor symmetries in this approach. 
For Majorana neutrinos G v has to be a direct product 
of two Z 2 , i.e. Z 2 < 8 » Z 2 to fully determine the mixing 
in neutrino sectoral]. For quarks, we have the same 
framework but since they are Dirac fermions, we do not 
have to be limit ed to Z 2 <g> Z 2 . In this case note that S T 
and C/J in Eqs. (6|7l should be replaced with 5^ and U^. 


II. A THEOREM FOR \U\ = \U\ T 

In this paper, as mentioned above, we define symmetric 
mixing matrix as 


\U\ = \U\ T (11) 

rather than the original definition of U = U T or the 
Hermitian case U = U'. The phenomenology is the same 
for \U\ = \U\ T and U = U T , but more general than the 
Hermitian case. 

Before we formulate our theorem that links geometrical 
properties of the flavor symmetry group to a symmetric 
mixing matrix, we will first define the geometric concepts 
which will be used. 

The Z 2 symmetries used in neutrino sector are actually 
just reflections or 180° rotations (the difference between 
them is trivial, any 180° rotation in 3-dimensional space 
can be changed to a reflection if we add an overall minus 
sign and vice versa). In 3-dimensional flavor space, 
going w.l.o.g. in the diagonal neutrino basis, the Z 2 
transformations correspond to putting minus signs to 
neutrino mass eigenstates: Vi -A — v t . The combined 
Z 2 (gi Z 2 in 3-dimensional flavor space corresponds to 
two reflections with respect to the direction of neutrino 
mass eigenstates. Or, in a picture we are more familiar 
with, if there are planes whose normal vectors are 
neutrino mass eigenstates, the Z 2 symmetries are just 
the mirror symmetries of those planes. Since Z 2 ® Z 2 
contains commutative mirror transformations or since 
the mass eigenstates are orthogonal, the mirrors should 
be perpendicular to each other, as shown in Fig. 0 by 
translucent squares. If the transformation from Ge is 
real (we will discuss the complex case later) in flavor 
space, it is an 50(3) transformation which can always 
be represented by a rotation. If the rotation axis is on 
the bisecting plane, which is defined as the plane that 
bisects the two mirror squares, or as the boundaries of 
octants in the diagonal neutrino basis, then we define 
that Ge bisects the two Z 2 ® Z 2 . We show two bisecting 
planes in Fig. |T] while all bisecting planes are shown in 



3 



¥ 


Figure 1. The geometrical relation of the mirror planes and 
the bisecting planes. The mirrors are placed on y — z, z — x 
or x — y planes. The two round disks are called bisecting 
planes because they bisect all the square mirror planes. The 
bisecting planes are boundaries of octants. 


Fig. [2] This gives now all definitions necessary for our 
theorem. 

Theorem A: 

If an 50(3) subgroup G contains two non-commutative 
Abelian subgroups G u and Gz, and if G u is isomorphic 
to Z 2 <S> Z 2 while Ge bisects the 7*2 <8> 72 , then G as a 
flavor symmetry can produce a symmetric mixing matrix. 

The definition of bisection and symmetric mixing are 
given previously. The subgroups G„ and Ge are required 
to be Abelian because the residual flavor symmetries are 
always Abelian |21U25] , and non-commutative so that the 
mixing is non-trivial. 

The proof of this theorem will be obvious after we in¬ 
troduce the general 50(3) rotation and the diagonaliza- 
tion below (see Eqs. ( pO| and (p4|). The Z 2 symmetries 
are applied to Majorana neutrinos and the bisecting rota¬ 
tion to charged leptons. However, one can also apply the 
theorem to quarks and obtain a symmetric CKM mixing. 
In the case of Dirac fermions, Z 2 is not necessary but suf¬ 
ficient. We will comment further on CKM mixing later. 
Because the axis of a bisecting rotation can be rotated on 
its bisecting planes, there are infinite bisecting rotations. 
Hence, Theorem A can produce infinite symmetric mix¬ 
ing matrices with one degree of freedom. Note that the 
unitary matrix with the constraint \U\ = \U\ T has only 
one prediction, see Eq. <§• 

Actually a lot of discrete flavor symmetries satisfy the 
conditions required by Theorem A, for example the tetra¬ 



Figure 2. The complete collection of all six possible bisecting 
planes and their geometrical relation with the mirror planes. 
A rotational symmetry with its axis on one of these bisecting 
planes can give, according to Theorem A, a symmetric mixing 
\U\ = \U\ T . 


hedral group T and octahedral group O which are just 
the widely used A 4 and 54 flavor symmetry groups, re¬ 
spectively (for geometrical interpretations on A 4 and 54 
see, for e.g. [25])• We can see in Fig. [ 3 ] that if we choose 
the three 180° rotational axes as x, y and 2 axes, which 
penetrate the tetrahedron through the two central points 
of two edges, then the bisecting planes are determined, 
as shown by dark blues circles. The tetrahedron is also 
invariant under the 120° rotation marked in Fig. [3] which 
is a bisecting rotation since the axis is on three bisecting 
planes. In explicit formulae, we say the tetrahedron is 
invariant under the following rotations, 

/o -1 0 \ 

Rbs = 0 0 -1 ; ( 12 ) 

\1 0 0 ) 


Si = diag(l, -1, -1); 5 2 = diag(-l, 1, -1). (13) 


Here Rb s is the 120° bisecting rotation and Si(i = 1,2) 
are the 180° rotations around the x and y axes. Eq. 
(121 can be obtained by requiring that l?b s (l,0,0) T = 
( 070 , 1 ) T , which means Rb s rotates the x axis to the z 
axis, as well as the other two relations l?b s (0,1,0) T = 
(-1, 0,0) T and R hs ( 0,0,1) T = (0, -1, 0) T . 

If Rbs and Si(i = 1,2) are the residual symmetries 
of the charged lepton sector and neutrino sector respec¬ 
tively, i.e. 


Rl s M e R hs = AI t ; SjM V S % = M v , (14) 

then according to Eqs. M > we can compute U( and U v 
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Figure 3. Tetrahedron symmetry. The dark blue circles show 
the bisecting planes. The axes of the 120° rotational symme¬ 
tries of the tetrahedron are on those planes, therefore accord¬ 
ing to Theorem A the tetrahedral group as a flavor symmetry 
can produce a symmetric mixing matrix \U\ = \U\ T . 


Figure 4. Octahedron symmetry. Similar to Fig. [3] accord¬ 
ing to Theorem A the octahedral group can also be used to 
produce a symmetric mixing matrix \U\ = |t/| T . 




from i?bs and «S). The result is 

, / 1 to u> 2 \ 

Ut = -j= I -1 -a; 2 -u> J ; U v = 1. (15) 

We see that U( is the Wolfenstein matrix, up to trivial 
signs. Therefore in this case the PMNS matrix U = UjU u 
is symmetric, i.e. \U\ = \U\ T . 

As another example we show in Fig. [4] that the 
octahedral symmetry, which is isomorphic to the widely 
used S 4 symmetry, also has the required properties for 
Theorem A. Figs. [3] and [4] show that the properties 
required by Theorem A are quite common in discrete 
groups with 3-dimensional irreducible real representa¬ 
tions. 

Now we present the theorem in explicit formulae. For 
simplicity, we choose a basis under which the mirror has 
a normal vector (1,0, 0) T , (0,1,0) T or (0,0,1) T so the 
mirror symmetry is just a reflection with respect to the 
y — z, z — x or x ~ y planes. Then the mirror transfor¬ 
mations through y — z, and z ~ x planes are 

S\ = diag(—1,1,1) and S 2 = diag(l, —1,1) 

respectively. Under this basis, the normal vectors n = 
(ni, n 2 , ns) T of the six bisecting planes satisfy one of the 
six conditions: 


\ n i\ = \rij |, (i,j = 1,2,3; i ± j ), (16) 

i.e. ni = ± 712 , n 2 = ±713 or n\ = ± 713 . The bisecting 
rotations i?bs with such an axis have special forms which 


we will show below. The neutrino and charged lepton 
mass matrices are now invariant under transformations 
of (Si, S 2 ) and i?bs respectively: 

SjM v Si = M v ; (i = l,2) (17) 

Rl s M e R hs = M t . (18) 

Diagonalizing these matrices with the transformation 


M„ = U V D„U? , M t = U 1 D 1 Uj 


(19) 


gives the PMNS mixing matrix U = lj)U v . According to 
Theorem A, U will be symmetric with proper ordering of 
the eigenvectors. 

Since we choose the basis in which Si and S 2 are di¬ 


agonal, M v is constrained to be diagonal by Eq. (171, 


hence C/„ is diagonal. As discussed at the end of the In¬ 
troduction, diagonalization of (Si,S 2 ) and i?bs will give 
Uy and Ui. So actually the key point of Theorem A can 
be stated as follows: 

For any SO( 3) matrix R, if R = i?bs there must be 
an unitary matrix U which is symmetric (\U\ = 
and can diagonalize R. The converse is also true which 
means R = Rb s is the necessary and sufficient condition 
for \U\ = \U\ T . 

Thus, in the basis we choose, we have a bisecting ro¬ 
tation to generate \U(\ = \Ui\ T and the mirror symme¬ 
tries to make Uy diagonal, therefore we get a symmetric 
PMNS matrix. In the above discussion we have explained 
the theorem in a specific basis, however, the physical re¬ 
sult is independent of any basis. One can choose another 
basis where the mirrors are not on the x — y, y — z and 
z — x planes, in which case the neutrino sector is not 
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diagonal and in general \Ui\ 7 / \U(\ T . However, the geo¬ 
metrical relation of the bisecting planes and the mirror 
planes makes sure that the product U}U V is symmetric. 


As for the explicit form of the bisecting rotation i?bs, 
we should first introduce the general rotation. The most 
general rotation in Euclidean space which rotates the 
whole space around an axis n = ( n\,n 2 ,n 3 ) T ( n.n = 1 ) 
by an angle <f> is 


R(n, <j>) 


nl+c (nl + n§) 
(1 - c)nin 2 - sn 3 
sn 2 + (1 — c)n\n 3 


(1 - c)n\n 2 + sn 3 
c + n 2 — cn 2 
—sn i + (1 — c)n 2 n 3 


-sn 2 + (1 - c)n 3 n 3 
sni + (1 — c)n 2 n 3 
c + n 2 — cn | 


where c = cos cf> and s = sin cf>. 


( 20 ) 


One can check that Eq. ( p0| does rotate the whole 
space around n by an angle while keeping n invariant. 
For example, when n = n z = (0,0,1) T we have 

c s 0 \ 

-scO, (21) 

0 0 1 / 


R(n z , <j>) = 


which is the familiar form of a rotation in the x — y plane 
around the z axis. 


For each of the six conditions in Eq. (16 1 we can get 
a bisecting rotation matrix from Eq. (20 1 . We use the 
symbol R^'^ to denote these bisecting rotations: 

R {±ij) = R{n\ ni=±nj ,</>). 

As an example, for n\ = n 3 , we have 


( d a p\ 

b h a , ( 22 ) 

q b d) 


where a = sni + (1 — c)nin 2 , b = (1 — c)n\n 2 — sni, 
d = c + nf(l — c). The remaining parameters p, q and h 
are determined by RR T = 1 if a, b, d are fixed. In general 
they are not equal to each other, but their precise forms 
are not important here. The point we should notice here 
is, if R has ni = n 3 , then the 12-element equals the 
23-element, the 21- the 32- and the 11- the 33-element. 


Conversely, if an SO( 3) matrix has the form of (22), then 


it must be a bisecting rotation with its axis on the x = z 
plane. This can be seen by solving (22 1 as an equation for 


(■ n,(f >) (the solution always exists since ( 20 ) contains all 
possible SO(3) matrices) and finding that the solutions 
always have ni = n 3 . 

R can be diagonalized by 


U^RUr = diag(e^,l,e-^) : 


(23) 


where the eigenvalues only depend on (j> while Ur only 
depends on n. As one can check numerically or by direct 
analytic calculation, \Ur\ has the following form, 


/ 1 — n\ 2 n\ 1 — n\ \ 

\U R \ 2 = -[ 1-nl 2n 2 2 l-n 2 2 ). (24) 

2 V 1 — 2 nl 1 -nj ) 


Here \Ur \ 2 is not but each element Xij of \Ur \ 2 is 

the absolute value squared of the ij-element of Ur. Note 
that the order of the columns in ( J24j ) can be changed since 
reordering of the columns of a diagonalization matrix is 
just a matter of permutation of eigenvectors. For n\ = 
n 3 , we recommend to write it in this order so that once 
one takes n\ = n 3 one immediately obtains a symmetric 
matrix. For the other cases such as n 3 = n 2 etc., we can 
always reorder the columns to get a symmetric matrix. 

From Eq. (|24l ), the proof of Theorem A is easy. One 
just sets Eq. p4j) equal to its transpose and finds nf = n 2 . 
There are other two possible permutation of the columns 
where (n\,n\,n\) T is the first or the last column of Ur, 
from which we can get = n 2 or n 2 = nf- 

This completes our proof of Theorem A. 

III. GENERALIZATION TO THE COMPLEX 
CASE 

The previous theorem only applies for flavor sym¬ 
metries with real representations, while some groups 
used in flavor symmetry model building enjoy complex 
representations. For the complex case, we cannot find a 
clear geometrical picture as was possible for real repre¬ 
sentations in 3-dimensional Euclidean space. However 
we can somewhat generalize the previous theorem to the 
complex case by finding some connections between the 
real and complex case |2B . In the following discussion, 
all unitary matrices are elements of SU (3) since the 
difference between 1/(3) and SU(3) is a trivial phase. 

Theorem B: 

If an SU(3) matrix T can be rephased to a real matrix 
R as follows, 

T = diag(e*“ 1 , e l “ 2 , e‘° 3 )R diag(e i/31 , e i02 , e i/3s ), (25) 

and if the R is one of the bisecting rotations with ni = nj 
from Theorem A, and if further 

&k +/3fc = 0(k ± i,j ), 

then T gives a symmetric mixing matrix 


( 26 ) 
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Note on Theorem B: k in Eq. (26) is the remain¬ 
ing number among {1,2,3} when \rii\ = \rij\ picks out 
two numbers for i and j. Since the rephasing matrices 
diag(e* Ql , e*“ 2 , e*“ 3 ) and diag(e 8/?1 , e l/3a , e*^ 3 ) should be in 
SU{3), it must hold oq + ct 2 + <^3 = Pi + P 2 + @3 = 0. So 
actually Eq. (26) is equivalent to at + a.j = — fit — Pj . 

As an exa mpl e, consider that the bisecting rotation is 
R( 13 ) in Eq. (221, then we have a 2 + P 2 = 0. In this case 
T is 


f f + ig ar]i pr ] 5 \ 

T (13) = bp 3 h aij 2 , (27) 

\ qvg b V4 f -ig ) 


could serve as a starting point or zeroth-order approxi¬ 
mation. 

The mixing scheme can be produced in A(96) group 
which can be defined by three generators a, b and c with 
the following properties m- 

a 3 = b 2 = ( ab ) 2 = c 4 = 1, 
caca = a~ 1 c~ 1 a = bcb ~ 1 = 6 _1 c&, 
cbc^b - 1 = bc- 1 b~ 1 c. (34) 

In a 3-dimensional faithful representation, a, b and c can 
be represented by [32] 


where rp are some phases, i.e. rp = 1. The 22-element 
is still h the (real) and 11-element is the conjugate of the 
33-element, as a result of a 2 + P 2 = 0. We also have 
r h r l' 2 V;iV 4 = 1 because oq + a 3 + Pi + /? 3 = 0. T^ 13 ) can 
be diagonalized by a unitary matrix which we call Uti 3 
and one can check that 


\Uri3 \ 2 — 


1 j— iu 1 — t 

2 ^ 2 Sip 1 

t __ g \ 

2 2 a v \ 

I 

1 f h—c v 

1 1 -c v 

1 -t 

, (28) 

1 9 1 t 

i + 9 J 

I 

2 2 s v 1 b 

2 ^2s v / 


= sin p and 

: ( — 1 + 2/ + h)/ 2 , 

(29) 

1 — h 

1 2(1 — Ccp ) 


(30) 


We can see that indeed |t/r 13 | = |f7n 3 | 


/o 1 o\ / 0 0-l\ 

a (3) =001, b (3) =0-10, 

\l 0 0/ \ -1 0 0/ 

/ i 0 0 \ 

c (3) = 0 -i 0 . (35) 

V° 0 

The mixing scheme is produced if a Z 3 subgroup gener¬ 
ated by T = a 2 (cb) 2 c and a Z 2 ® Z 2 subgroup generated 
by S\ = b, S 2 = a 2 c 2 ac 2 is applied to charged leptons 
and neutrinos respectively. To be precise, 

/ 0 0 i \ 

Z 3 : T = -1 0 0 (36) 

\ 0 ! 0 / 

and the two Z 2 are generated by 


IV. APPLICATION 

In this section we will apply our theorems to an ac¬ 
tual mixing scheme. After the T2K neutrino experiment 
measured a large non-zero 9 13 in 2011 [29], many models 
have been proposed to explain the result. Refs. Esum 
scanned a series of discrete groups (A(6n 2 ) and T^v), and 
one of the found schemes was quite close to the T2K re¬ 
sult at the time. In the standard parametrization, the 
angles are 

023 = 012 = tan" 1 -J— = 36.21° , (31) 

v3 +1 

and 

(32 > 

In total the PMNS matrix is 

/H3 + V3) Gs 

u = -73 73 71 • ( 33 ) 

U( 3 "G3) |( 3 + V3); 

While this mixing scheme is ruled out by current data, it 
fulfills our criterion of a symmetric mixing matrix, and 


/o 0 l\ /-1 0 0 \ 

Z 2 ® Z 2 : Si = - 0 1 0 , S 2 = 0 1 0 . 

\ 1 0 0 / \ 0 0 -1 / 

(37) 

From our theorem, it is easy to see that this can produce 
symmetric mixing. In the diagonal neutrino basis we 
transform T to T d , 

( i _l_ i \ 

2 v /2 2 \ 

_ 7f 0 72 (38) 

i _j_ _ i I 

2 72 2 / 

and (Si, £ 2 ) to (Si d , S 2d ): 

Sid = diag(—1, -1,1); S 2d = diag(-l, 1, -1). (39) 

Then, according to our theorem, we see that T d 
can be rephased via a transformation defined in Eq. 
p5} to a bisecting rotation R with a 2 = P 2 = 0, 
(a7,a 3 ) = (7r/2, —7 t/ 2) and {pi,P 3 ) = (0,0). Or, in a 
simpler way, T d has the form of Eq. ( |27[ ) . So the mixing 
matrix should be symmetric. 

A dynamical realization of the mixing scheme in A (96) 
has been studied in Ref. [33]. That model is rather com¬ 
plicated using both 3- and 6-dimensional representations. 
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Here we present a simpler model which only uses two ad¬ 
ditional sets of scalar fields ft, ft, and features all par¬ 
ticles in the same 3-dimensional representation of A(96) 


of (361 and (37). 


t,£ c ,v , ft, ft ~ 3. 


(40) 


We use the representation in (361 and 


rather than 


(381 and (391 because the Clebsch-Gordan (CG) coeffi¬ 
cients are simpler. The result does not depend on the 
basis. The CG coefficients we will use in this representa¬ 
tion are 


3 0 3 1 : (41) 

3 0 3 0 3 ->• 1 : e ijk (42) 

3030303—^1: S ijmn (43) 

3030303—>• 1 : S im 8 jn S in S jm . (44) 


Here etj k is the Levi-Civita tensor (or order 3 antisym¬ 
metric tensor) and ftj mn is defined as 


(i = j = m = n), 
otherwise. 


(45) 


The invariant Lagrangian in the lepton sector is 

C = y[e ijk 4>%£t + V 2 8 ijmn ft ft £m£-n + 

ft Un T V'2 8rmVjn ( Pi ft Vn ■ (46) 

After symmetry breaking, ft and ft obtain the following 
VEVs: 

(ft) = ^(1,-1, i) , (ft) = V*(l, 0,1). (47) 

In the charged lepton sector Aft = is 

( u x + iy —ix — y \ 
x — iy u ix — y I , (48) 

ix — y —ix — y u J 

while the neutrino mass matrix is 

( A 0 B\ 

M v = ft| 2 0 0 0, (49) 

\ B 0 A j 

where A = ft + ft; B = ft with 

u = 2|ftl“ + 21 2 /| | “ + \y l 2 + 2/11 2 (50) 


x = \y[\ 2 -3\yi\ 2 -2Re[yi*yi\, (51) 


y = 2R e[y{*y e 2 \. (52) 

Mf and M v can be diagonalized by the following unitary 
matrices 


U v 


J_ 0 A 

02 U 02 

0 10 

J_ o -A- 
02 U 02 


(53) 


1 


1 


U * = - 1 -u 2 

v ^ \ —i —i —i 


(54) 


Thus, the PMNS matrix is 


2+Cc! 

cA 

— i-\-uj 

06. 

03 

06. 

1+4 

1 

l — l 

V6 

03 

Vs 

l+UJ- 

L0 


06 

03 

06 


It is related to the matrix in Eq. (j33J) via 

^pmns = diag(e I/3l ,e ,/32 ,e I/33 )t/diag(l,e z 


(55) 


(56) 


where ft = 105°,ft = 225°, ft = 165°, = 45°, 

«2 = 90°. Here a\, a 2 would be the Majorana phases if 
the couplings y\ and ft i n Eq. (46) were real. Therefore 
even though the Dirac-type CP is conserved in this 
model, generally there is still CP violation due to the 
non-zero Majorana phases, unless the phases of ft and 
ft are tuned to exactly cancel a±, « 2 - 


Our theorem can also be applied to the quark sector. 
One just assigns the bisecting rotational symmetry to the 
residual symmetry of up-type quarks (or down-type) and 
the mirror symmetries to that of down-type quarks (or 
up-type), then the CKM mixing will be symmetric. How¬ 
ever building a realistic model for the CKM mixing is a 
somewhat more difficult task. Compared to the lepton 
sector where hundreds of flavor symmetry models have 
been proposed, for the quark sector much fewer models 
exist. This is due to the fact that the small CKM mixing 
angles do not have straightforward geometric interpreta¬ 
tion, which is the basis of discrete flavor symmetry build¬ 
ing. Among the existing models for the CKM mixing, we 
cannot find one that fulfills our criteria (exceptions are 
of course the trivial cases in which one interprets the 
CKM matrix as the unit matrix or as a matrix which 
only consists of the Cabibbo angle), and scanning all dis¬ 
crete groups for the flavor symmetry of quarks is out of 
the main purpose of this paper. Anyway, when looking 
for flavor groups to build models for the quark sector, 
our theorem could be a guidance because when a mix¬ 
ing scheme generated from a flavor symmetry is close to 
realistic CKM mixing, then it must be also close to a 
symmetric form. 


V. CONCLUSION 

A possible zeroth-order, but surely aesthetically at¬ 
tractive, mixing Ansatz for the CKM and PMNS matri¬ 
ces is that they are symmetric. The origin of symmetric 
PMNS and CKM matrices from the viewpoint of flavor 
symmetry models has been the focus of our paper. 

We have proposed a theorem on the relation between 
symmetric mixing matrices and geometric properties of 
discrete flavor symmetry groups. An illustrative connec¬ 
tion between the rotation axes of the geometric body 











associated to the symmetry group exists, and shows that 
popular subgroups of SO( 3) such as A 4 and 64 can lead 
to symmetric mixing matrices. 

Groups with complex irreducible representations do 
not easily allow for a geometrical interpretation, but a 
partial generalization of our theorem is possible, which 
can then apply to SU(3) subgroups such as A(96). A 
previously studied mixing scheme that turns out to cor¬ 
respond to a symmetric PMNS matrix was used as an 
explicit example. 

The connection of geometric properties of discrete 


groups and possible features of the mixing matrices may 
have further applications. 
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